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I.  INTRODUCTION 


The  determination  of  the  charged  particle  distribution  in  the  plasma 
sheath  which  envelopes  a  reentry  vehicle  is  important  in  understanding  the 
electromagnetic  response  of  the  vehicle.  Hence,  the  Trailhlazer  Reentry 
Test  Program  was  established  at  AFCRL  to  allow  direct  measurements  of  the 
electron  density  distribution  in  the  plasma  sheath. 

This  report  is  an  analysis  of  ion- collecting  electric  probe  data  ob¬ 
tained  between  altitudes  of  220  Kft  and  270  Kft  for  the  third  flight  of 
Trailblazer  II.  The  blunt-nosed  vehicle  had  a  nose  radius  of  .527  ft  and  a 
9°  conical  afterbody.  Its  angle  of  attack  during  reentry  was  approximately 
10°  and  it  was  spinning  about  its  center  line  axis  at  approximately  10  revo¬ 
lutions  per  second.  The  ion  collection  data  were  obtained  using  three  £- inch 
diameter  flush-mounted  copper  electrodes  that  were  insulated  from  the  rest 
of  the  vehicle  surface.  The  probes  were  biased  at  voltages  ranging  from 
-5  to  -30  volts  and  were  located  at  the  nose  cap  center  point  and  at  S/R  = 
.^75  and  2.85  respectively.  Here,  S  is  arc  distance  back  from  the  nose  cap 
center  point  and  R  is  the  nose  radius.  Positive  current  traces  for  these 
three  probes  are  shown  in  Figs.  1  through  4. 

In  analyzing  probe  data,  one  can  either  solve  the  direct  problem, 
whereby  one  assur.es  the  flow  and  chemistry  to  be  known  and  then  predicts 
-.he  response  of  the  probe  to  a  given  applied  voltage,  or  the  indirect  prob¬ 
lem,  where  one  tries  to  infer  the  chemistry  frctn  the  probe  current-voltage 
characteristics.  Each  method  has  its  advantages  and  disadvantages.  With 
the  direct  approach,  one  may  begin  with  an  incorrect  ion  number  density  pro¬ 
file  and  therefore  not  be  able  to  predict  the  probe  response  even  if  the 
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mechanism  by  which  ions  are  transported  to  the  probe  is  correctly  formulated. 
With  the  indirect  approach,  there  is  the  problem  of  associating  a  unique 
chemistry  model  with  the  given  probe  characteristic.  In  either  case,  there 
is  the  possibility  that  either  the  chemistry,  the  modeling  of  the  ion  sheath 
structure,  or  both  may  be  incorrect.  In  this  study,  the  so-called  direct 
approach  will  be  utilized. 

To  obtain  a  theoretical  prediction  of  the  current  as  a  function  of 
altitude,  body  position,  and  applied  potential,  the  charged  particle  dis¬ 
tribution  in  the  boundary  layer  must  be  known.  To  this  end,  a  laminar  clean 

air  chemistry  model,  in  which  either  vibrational  equilibrium  of  vibrational 

1  2 

nonequilibrium  could  be  assumed,  was  used  ’  .  Although  it  is  possible  that 
discrepancies  between  theory  and  experiment  may  be  caused  by  lack  of  knowledge 
about  the  clean  air  chemistry,  the  approach  taken  in  this  study  is  to  initi¬ 
ally  assume  the  clean  air  chemistry  is  known. 

Theories  describing  ion  collection  by  electric  probes  can  be  very  com¬ 
plex  depending  on  the  values  of  such  parameters  as  the  ratios  of  sheath 
thickness  £ _  to  ion-neutral  mean  free  path  A^_n  or  boundary  layer  thick¬ 
ness  V,  „  .  If  the  former  is  large,  the  effects  of  collisions  inside  the 

b.l. 

sheath  must  be  taken  into  consideration.  If  the  latter  is  not  much  less 
than  one,  effects  of  convection  within  the  sheath  may  have  to  be  considered. 

To  simplify  the  data  analysis,  the  experiment  was  designed  with  the 

probe  voltages  sufficiently  large  in  magnitude  that  the  full  random  flux  of 

ionized  particles  of  one  sign  would  be  collected  at  the  sheath  edge.  Then 
1  1 

a  theory  ’  first  proposed  by  Bredfeldt  and  Scharfman  may  be  applied  to 

■3  q 

analyze  the  data.  This  theoryJ  was  used  by  Hayes  to  infer  electron  number 
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densities  at  a  number  of  altitude  and  body  station  combinations  for  flight 
number  three.  Except  at  an  altitude  of  27O  Kft,  there  was  good  agreement 
between  electron  number  density  inferred  from  the  data  and  that  predicted  by' 
flow  field  calculations  in  the  nose  cap  region  of  the  vehicle.  However,  the 
agreement  on  the  conical  portion  of  the  vehicle  was  not  good. 

Because  of  this  and  the  dichotomy  in  the  Br edf eld - Scharfman  theory  of 
simultaneously  assuming  a  collisionless  flux  of  particles  into  the  outer 
edge  of  the  sheath  and  a  collision-dominated  sheath,  it  was  decided  to 
examine  other  approaches.  As  explained  in  detail  in  this  report,  these 
include  a  load-line  technique  that  accounts  for  collisions  in  the  sheath  and 
the  large  changes  in  density  that  may  occur  across  the  sheath  and  an  analy¬ 
sis  of  convection  as  a  possible  ion  transport  mechanism  within  the  sheath. 
The  possibility  of  current  enhancement  because  of  charge  depletion  near  the 
outer  edge  of  the  sheath  is  also  examined. 

If  the  applied  potential  is  sufficiently  high  that  the  sheath  is  many 
mean- free  paths  in  thickness ,  the  ion  drift  velocity  may  be  governed  by  a 
mobility-limited  drift  relationship^.  But,  this  relationship  in  only  valid 
if  the  energy  an  ion  gains  in  one  mean  free  path  is  much  less  than  its  ther¬ 
mal  energy.  Keeping  this  limitation  in  mind,  it  was  decided  to  use  these 
other  approac les  in  conjunction  with  a  mobility- limited  ion  drift  to  try  to 
achieve  bette  .•  agreement  between  the  data  and  flow  field  predictions  on  the 
conical  portion  of  the  vehicle  at  all  altitudes  considered  and  on  the  nose 
cap  at  270  Kft. 
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II.  EQUATIONS 


If  the  ion  sheath  is  thin  with  respect  to  the  boundary  layer,  convec¬ 
tion  in  the  sheath  may  he  neglected  and  a  one-dimens ionai  analysis  may  he 
used  .  However,  if  the  sheath  is  sufficiently  thick  that  convection  can't 
he  neglected  within  it,  two-dimensional  effects  must  he  included J’  .  Since 
the  prohe  diameter  is  slways  considerably  smaller  than  the  boundary  layer 
thickness,  if  convection  is  important,  fringing  effects  will  also  he  signi¬ 
ficant.  Hence,  the  problem  (cf.  Fig.  5)  is  not  only  x  dependent  and  y  de¬ 
pendent  hut  also  z  dependent,  clearly  a  most  complicated  situation. 

Even  if  one  assumes  that  the  z  dependence  can  he  approximately  accounted 
for  by  applying  a  correction  factor  to  a  solution  that  assumes  g  .  =  0,  the 

■V 

problem  is  still  forbidding.  The  full  equations  describing  the  motion  of  the 
electrons  and  ions  over  the  electrode  are,  assuming  a  steady  state, 

V\ 

Ion  Cont.  V'N.  (U  +  K ,  X7  V  -  D,  - )  =  0  (l) 

l  i  y  1  w- 


Electron  Cont.  V,We  (u  "  KeV  v  -  De 


(2) 


Poisson 


(3) 


vhere 


V  is  the  potential 

is  the  ion  number  density 
Ne  is  the  electron  number  density 
f 0  is  the  dielectric  constant  for  a  vacuum 
U  is  the  velocity  of  the  neutral  particles 
Ik  is  the  ion  diffusion  coefficient 

De  is  the  electron  diffusion  coefficient 

Kk  is  the  ionic  mobility 

Ke  is  the  electron  mobility 

It  is  assumed  here  that  the  gas  is  only  slightly  ionized  so  the  velo¬ 
city  and  temperature  profiles  foi' the  neutrals  can  be  determined  independently 
of  the  species  conservation  equations. 

Inside  the  ion  sheath,  Ng  --  0  en 

N.  =  1-^-2  V2V 
1  e 

r  vv> 

(U  +  %  VV  -  Di  -  0 

V*  (U  +  K.^V  -  D.  ?  ^  =  0  00 

Equation  (4)  is  a  fourth  order  partial  nonlinear  equation  in  the  poten- 
tiex  V.  Cumbersome  numerical  techniques^  have  been  employed  in  an  attempt 
to  intergrate  Eqn.  (4).  Dukowicz simplified  Eqn.  (4)  somewhat  by  neglect¬ 
ing  the  underlined  diffusion  term.  He  solved  the  problem,  numerically  for  an 
incompressible  inviscid  flow  past  an  ion- collecting  flat- plate  probe.  His 
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solutions  are  valid  for  a  constant  convective  velocity  and  ion  number  density 
upstream  of  the  electrode.  Hence,  they  are  not  applicable  to  a  ■'"Cw.itry  boun¬ 
dary  layer  problem. 

In  the  regime  where  convective  and  charge  depletion  effects  are  small, 
but  where  property  variations  across  the  sheath  are  significant,  one  can  uti- 

Q 

lize  a  so-called  "load-line"  technique  to  predict  the  response  of  the  probe 

to  a  given  applied  potential.  The  load  line  equation  is  derived  from  Eqn. 

(*0  by  letting  —5? .  =  — ~  =  0,  and  neglecting  convection  and  diffusion. 

z 

Then  the  flow  of  ions  is  one-dimensional  with  the  ion  continuity  equation 
becoming  (cf.  Fig.  6) 

j=eH.v.  (5) 

where  the  velocity  of  the  ions  v^  =  K^E  is  toward  the  wall  as  is  the  ion 
current  density  j . 

Combining  Eqn.  (5)  with  Poissons  Equation: 


N, 


l 


So  dj 

e  dy 


results  in  the  differential  form  of  the  load  line  equation 


(6) 


With  j  e=  JTT  being  constant  across  the  sheath.  Integrating  Eqn.  (j)  across 
LL 

•che  sheath  and  applying  the  bound1  ry  conditions  V(  fs)  =E  (  fs)  =0,  V(0)  =  V  , 
leads  to  the  load-line  equation: 
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where  E  Lg  the  electric  1'ield.  Since  j  is  constant  for  a  r‘  ,hn  Eqn.  (6) 

can  he  solved  for  ,j  in  terns  of  a  double  integr*’,’  ever  the  sheath.  For  a  given 
Vw  and  it  describes  the  ability  u.n  electric  sheath  to  transport  a  cer¬ 

tain  flux  of  ions  by  means  of  a  mobility- limited  drift.  By  maintaining  the 
option  of  allowing  u>  vary,  the  large  variation  in  properties  through  the 
boundary  layer  are  taken  into  account. 

A  Fortran  listing  for  the  computer  code  that  both  computes  the  load-line 
current  J  and  other  relevant  quantities  may  be  found  in  the  appendix. 

At  a  given  altitude  and  body-station  combination,  the  following  profiles 
from  a  boundary  layer  program  are  inp.it  as  functions  of  distance  y  out  from 
the  wall.  They  are  density  p(RHO),  temperature  ratio  T/Te  (THETA.),  ion  = 
electron  number  density  (ED),  and  velocity  ratio  u/uc  (FP).  Here,  T  refers 
to  temperature,  u  to  velocity  in  the  x  direction,  and  the  subscript  "e"  refers 
to  conditions  at  the  outer  edge  of  the  boundary  layer.  The  particular  boundary 

p 

layer  profiles  used  for  this  study  were  obtained  by  Lew  of  GE.  His  profile 
data  are  given  at  thirty  unequally  spaced  points  throughout  the  boundary  layer. 
To  obtain  'alues  of  this  data  at  other  y  coordinates  in  the  boundary  layer,  sub¬ 
routine  GETRH0  is  used.  In  this  subroutine,  a  Lagrangian  three-point  interpola¬ 
tion  scheme  is  used  to  find  values  for  the  dependent  variable  at  the  desired 
mesh  points. 

The  convective  current  density  jc  is  found  simply  from  the  product 

jc  =  e  u  N.  (9) 

The  integrated  convective  current  density  is  found  from 

ye 

Jc  =  ^  £  Jc  dY  (io) 
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where  dp  is  the  prohe  diameter  «=  .63  cm.  The  integration  is  performed  by 
using  Simpson’s  rule. 


The  diffusion  current  density  Jg  is  obtained  from 


=  e  D, 


i  dy 


dN 


where  the  derivative  — is  found  by  using  a  central  difference  formula 
and  the  diffusion  coefficient  is  found  from  the  Einstein  relation 


Ki  '  It 


and  a  curve  fit  for  Mason's  calculations  for  the  mobility 


%  - 11  A/jP 


v-sec 


The  constant  density  load-line  current  is  found  by  first  of 

— —  —1 

all  finding  the  average  density  (y)  dy  throughout  the 

sheath,  using  Eqn.  (v)  to  calculate  the  average  mobility  and  then  using 
Eqn.  (9)  to  find 

The  bulk  of  the  main  program  is  used  to  find  the  load-line  current  den¬ 
sity  JLL*  To  do  this,  Eqns.  (k),  (6)  and  (13)  are  written  in  the  form 


A- 

dy 


-  dV 


For  a  given  sheath  thickness  £s  .  These  equations  are  integrated 

in  toward  the  wall  using  E  (Q  =v  (4)  =  0  as  initial  conditions.  To  per¬ 
form  this  integration,  the  subroutine  MERKUT  is  used  which  is  Merson's  adapta¬ 
tion  of  a  Runge-Kutta  integration  scheme.  It  was  found  convenient  to  work  in 
terms  of  the  variables  E  2  and  V  2/3  where  (  )  denotes  division  by  4 2  J  ,  in 
performing  this  integration.  With  these  variables,  JLL  did  not  appear  expli¬ 
citly  in  the  equations.  Upon  reaching  the  wall,  is  found  from 


JlL  2  V2  (0) 


and  E  and  V  are  found  from 


V  = 

The  percentage  error  (as  computed  by  Merson’s  method  )  in  E2  and  j_s  held 

to  less  in  one  part  in  10^.  B;y  picking  a  number  of  different  values  for  $ 

s 

expending  from  0  to  y  ,  a  complete  load-line  curve  can  be  computed. 

D  •  X  « 

To  determine  the  importance  of  including  the  proper  boundary  layer  den¬ 
sity  profile  in  the  sheath  equations,  a  second  load-line  calculation  was  per¬ 
formed  with  the  density  kept  constant  at  its  average  value  throughout  the 
sheath.  For  this  special  case  of  K±  constant,  Eqn.  (8)  reduces  to  the  contin¬ 
uum  equivalent  of  the  Child- Langmuir  law 

JLLCD  =  l  £0  Ki  “7T  (l6) 


where  J  has  units  of  current  per  unit  area, 
LLCD 


If  convective  effects  are  important  and  diffusion  is  negligible,  the 
Quasi-One-Dimensional  (Q-l-D  )  approach  discussed  in  Ref.  9  may  be  used  to 
provide  a  rough  estimate  of  the  current  collected  by  the  probe.  The  basic 
idea  here  is  that  the  ions  are  convected  into  the  sheath  (cf.  Fig.  6)  and 
then  attracted  to  the  vail  by  the  mobility-limited  drift  phenomenon.  The 
theory  is  only  valid  when  flow  quantities  vary  much  more  strongly  in  the  y 
*.han  in  the  x-direction.  But,  we  will  apply  it  to  the  problem  at  hand  anyway 
as  a  means  of  estimating  convective  effects. 

We  first  use  Eqn.  (8)  to  find  the  load-line  current  density.  The  aver¬ 
age  convective  current  density-  per  unit  probe  area  Jc  is  found  from 


0 


where  dP  is  the  electrode  diameter.  Setting  J  as  calculated  by  Eqns.  (8) 
and  (17)  equal  to  one  another  allows  us  to  calculate  a  Q-l-D  current  density 
•Jq-I-P  and  sheath  thickness. 
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IH.  ANALYSIS  OF  DATA 


The  load-line  program  has  been  used  to  calculate  the  load-line  current 
Jpp,  constant  density  load  line  current  JTT^n,  diffusive  current  J^,  convec¬ 
tive  current  density  «jc,  and  integrated  convective  current  JQ  at  ten  differ¬ 
ent  altitude-body  station  combinations.  Results  are  shown  on  Figs.  7-19. 

The  label  at  the  top  of  each  figure  indicates  the  altitude  and  body  station 
combination  of  interest.  In  each  case  the  integrated  convective  current  Jc 
attains  its  maximum  value,  of  course,  at  the  outer  edge  of  the  boundary  layer. 
This  asymptotic  value  is  the  maximum  current  the  probe  is  capable  of  collect¬ 
ing  a^s’mung  fringing  effects  are  unimportant  in  calculating  the  convective 
current  collected.  The  convective  current  density  jc  attains  a  peak  -value 
within  the  boundary  layer  because  the  ion  number  density  approaches 
zero  near  the  outer  edge  of  the  boundary  layer  and  the  convective  velocity 
u  approaches  zero  near  the  wall. 

At  all  altitudes  the  diffusion  current  neglecting  charge  depletion 

effects  J,  is  at  least  an  order  of  maimitude  less  than  the  convective  eur- 
d  ~ 

rent  Jc  that  is  eonvected  through  the  entire  boundary  layer.  At  that  place 
in  the  boundary  layer  where  reaches  its  peak  value,  the  diffusion  current 
approaches  zero. 


The  load-line  current  density  or  increases  as  the  sheath 

thickness  decreases  in  accordance  with  Eqns.  (8)  and  (l6) .  Since  much  of 

the  experimental  data  was  obtained  at  probes  voltages  of  -5  and  -15  volts, 

theoretical  results  are  shown  at  an  intermediate  value  of  -10  volts.  The 

2 

results  can  be  scaled  for  other  voltages  simply  by  following  a  Vv  scaling 
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law.  Although  only  shoved  explicitly  on  Fig.  10,  is  generally  about 


30$  lower  than  JTT  over  the  regime  studied. 


Experimental  results  for  currents  collected  as  functions  of  altitude 

are  shown  on  Figs.  1-4  for  V  =  -15  volts.  Results  of  V  =  -5  volts  are 

w  w 

generally  around  30$  lover.  Except  for  the  probe  located  at  the  nose  cap 
center  point  (S/R  =  0) ,  all  shoved  roughly  a  30$  to  50$  fluctuation  with 
angle  of  attack.  Since  a  probe  has  a  surface  area  of  .316  cm  ,  the  values 

must  be  multiplied  by  3*l6  to  be  converted  into  current  densities  having 

.  o 

units  of  amps/cnf. 


Summary  plots  of  both  the  experimental  and  theoretically  predicted  cur¬ 
rent  densities  as  functions  of  arc  position  along  the  body  S  are  shown  in 
Figs.  7,  8  and  9.  The  experimental  results  are  obtained  from  Figs.  1-4 
with  the  assumption  that  the  peak  value  of  Jexp  in  a  given  (small)  altitude 
range  corresponds  to  the  probe  positioned  on  the  windward  side  of  the  vehicle. 
Conversely,  the  minimum  value  is  obtained  when  the  probe  is  on  the  leeward 
side  of  the  vehicle  at  the  maximum  angle  of  attack  <Y  =  10°.  Hence, 
the  nose  cap  center  point  location  corresponds  to  S  =  .09  ft,  the  probe 
located  at  S/R  =  .475  corresponds  to  S  =  .12  feet  (windward)  and  .34  feet 
(leeward) ;  and  the  probe  located  at  S/R  =2.58  correspond  to  S  =  I.25  feet 
(windward) and  1.47  feet  (leeward). 


Except  for  currents  measured  by  the  probes  located  on  the  spherical 
nose  (s/R  =  0  and  .475)  at  altitudes  close  to  270  Kft,  all  the  experimental 
currents  were  less  than  the  total  current  being  convected  through  the  boun¬ 
dary  layer.  Here  total  current  is  defined  as  that  being  convected  through 
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a  cross-section  perpendicular  to  the  probe  having  a  width  equal  to  the  probe 
diameter  and  a  height  equal  to  the  boundary  layer  thickness.  Since  Hie  boun¬ 
dary  layer  is  always  many  probe  diameters  in  thickness,  and  the  probe  is  only 
as  wide  as  one  probe  diameter  at  one  point,  this  should  definitely  be  an 
upper  bound  on  total  current  being  colle  cted  by  the  probe .  Even  if  the  ion 
sheath  were  to  extend  to  the  edge  of  the  boundary  layer,  not  all  the  ions 
entering  the  sheath  would  re  collected  because  of  the  way  the  electric  field 
and  convective  velocities  determine  the  ion  trajectories. 

By  comparing  Jc  with  Jexp,  one  can  conclude  that  the  predicted  ion  den¬ 
sity  profiles  are  plausible  except  between  S/R  =  0  and  S/R  =  .^75  at  z  =  2J0 
Kft. 


The  diffusion  currert  density  is  always  less  than  Jg  .  A  strictly 
one-dimensional  analysis  if  the  sheath  is  in  a  collision-dominated  regime 
would  involve  use  of  the  relation  to  compute  the  sheath  thickness 

and  hence  the  current  density.  This  results  in  relatively  thin  sheaths  and 
current  uensities  at  least  an  order  of  magnitude  less  than  those  measured. 
The  corresponding  curves  are  labeled  JLL-d  and  shown  on  Figs.  7>  8  and  9* 

At  270  Kft,  is  always  greater  than  J^.  In  other  words,  the  diffusion 
process  never  supplies  enough  current  to  satisfy  the  load-line  equation  if 
charge  depletion  effects  are  neglected.  Hence,  at  270  Kft  the  diffusion 
current  is  plotted  vs.  body  station  S  in  Fig.  8.  At  each  S,  the  plateau 
value  of  is  chosen. 

One  possibly  significant  effect  that  has  not  been  investigated  in  de¬ 
tail  is  the  phenomenon  of  local  charge  depletion  in  the  boundary  layer.  As 
ions  are  collected  by  the  probe  the  ion  number  density  immediately  outside 
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the  sheath  region  (and,  of  course,  the  electron  number  density)  will  he 
reduced  somewhat  as  ions  are  removed  from  the  boundary  layer.  This  local 
1 eduction  in  %  will  tend  to  increase  over  that  which  is  predicted  by  the 
boundary  layer  program. 

If  charge  depletion  is  important,  then  the  convective  and  diffusive 
terms  in  Eqn.  (l)  are  certainly  important  in  describing  flow  phenomenon 
in  the  charge  depletion  region.  Inside  the  sheath,  the  diffusion  term  may 
be  important  near  the  outer  edge  and  should  not  be  neglected  in  the  load- 
line  equation.  For  sheaths  sufficiently  thick  with  respect  to  the  boundary 
layei ,  the  convective  term  will  be  important  also  inside  the  sheath. 

To  estimate  the  order  of  this  effect,  the  width  2wd  of  this  diffusion 
region  may  be  estimated  from 


2“a  ■  2  /  VHr 


(18) 


where  D.  is  the  ion  diffusion  coefficient  at  the  outer  edge  of  the  sheath 
i 

and  t  is  the  flow  transit  time  over  the  probe.  The  rate  at  which  ions 
tr 

are  convected  into  this  region  is 


2  e  N.  u  w.  d 

i  d  p 


(19) 


where  d^  =  .635  cm  is  the  probe  diameter. 
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If  this  current  is  continually  being  attracted  to  the  probe,  the  net 


increase  in  probe  current  density  vill  be  approximately 


e  N.  u 


(20) 


For  body  stations  on  the  nose  and  on  the  cone  at  220  Kft  and  270  Kft, 
this  quantity  was  evaluated  and  found  to  be  about  3  on  the  nose  and 
10  on  the  cone.  The  width  of  the  diffusion  region  was  found  to  be 
of  the  order  of  d^.  At  220  Kft  there  is  the  possibility  of  the  neces¬ 
sary  ion  current  being  supplied  by  charge  depletion  effects.  At  270  Kft 
the  total  current  in  the  nose  region  is  still  too  low  to  match  the  data 
out  the  current  due  to  convection  and  charge  depletion  may  be  sufficient 
on  the  cone. 

More  detailed  calculations  need  to  be  performed  to  clarify  this 
point.  The  results  of  preliminary  calculations  discussed  here  are  shown 
as  solid  circles  on  Figs.  7  and  3. 

Although  fringing  effects  undoubtedly  would  not  increase  the  upper 
bound  previously  calculated  for  Jc,  the  effective  surface  area  of  the 
probe  can  be  assumed  to  be  enlarged  by  one  quarter  of  a  torus  having  a 
centroid  at  the  center  of  the  electrode.  The  correction  factor  by 
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which  the  actual  surface  area  of  the  electrode  most  "be  multiplied  "to 
give  an  effective  collection  area  for  diffusion  cuirent  is  then 


c 


1  + 


+  2  (■ 


<fs 


r 


P 


2 

) 


(21) 


This  equation  is  useful  only  if  one  can  calculate  a  sheath  thick¬ 
ness  $~s.  As  described  earlier,  in  a  region  where  convective  effects 
dominate  over  diffusive  effects  in  transporting  ions  into  the  sheath, 
the  Q-l-D  approach  may  "be  used  to  obtain  an  estimate  of  the  sheath 
thickness.  Results  for  current  density  obtained  by  this  method  shows 
disagreement  with  the  theory  by  an  order  of  magnitude  on  the  nose. 

Back  on  the  cone,  there  is  reasonable  agreement  between  theory  and 
experiment  at  270  Kft  while  the  Q-l-D  current  is  too  low  by  a  factor 
of  3  "to  7  at  220  Kft.  The  Q-l-D  results  plotted  here  are  for  Vw  =  -10 
volts.  If  we  had  used,  instead,  results  for  Vw  =  -15  volts,  the  pre¬ 
dicted  currents  would  have  been  twice  as  high  and  hence  closer  to  the 
experiment  current.  The  theoretically  predicted  voltage  variation 
does  not  agree  with  the  exp_rimentally  measured  current  increase 
of  roughly  a  factor  of  30$  as  Vw  went  from  -5  to  -15  volts.  Hence, 
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it  would  be  misleading  to  claim  agreement  between  theory  and  experiment  if 
we  could  match  the  data  at  one  (higher)  voltage.  This  is  why  a  representa¬ 
tive  voltage  of  10  volts  was  picked  for  this  analysis. 

Since  the  Q-l-D  method  obviously  is  not  always  applicable,  and  the  dif¬ 
fusion  current  without  charge  depletion  is  too  small  to  match  the  experimental 
current,  a  sheath  thickness  was  estimated  by  setting  Jc  =  JeX£.  The  results 
for  a  few  body-station  altitude  combinations  are  shown  in  Fig.  21.  They 
indicate  values  of  c  in  Eqn.  (21 )  from  about  3  to  30  with  the  larger  values 
being  associated  with  the  higher  altitudes.  If  convection  is  a  significant 
transport  mechanism,  these  results  indicate  that  fringing  effects  can  be 
quite  important  in  influencing  the  total  diffusion  current  collected.  Of 
course,  in  a  regime  where  fringing  effects  dominate,  the  utility  of  Eqn. 

(2l)  becomes  very  limited. 

It  should  be  noted  that  with  all  of  the  continuum  theories  discussed 
here,  the  sheath  is  many  ion-neutral  mean  free  paths  in  thickness.  With 
sheaths  of  these  thicknesses,  it  would  not  be  proper  to  use  the  thermal 
drift  idea  of  Ref.  3  to  calculate  the  ion  current  flowing  into  the  outer 
edge  of  the  sheath.  However,  in  Ref.  1,  the  calculated  sheath  thicknesses 
were  only  a  few  mean  free  paths  or  less  in  thickness.  Hence,  it  might  have 
been  consistent  in  that  case  to  use  a  thermal  drift  velocity  to  calculate  an 
ion  flux  into  the  sheath. 

The  main  inconsistancy  in  using  the  thermal  flux  plus  the  mobility- 
limited  drift  relationship  simultaneously  is  that  they  lave  different 
regions  of  the  validity.  The  former  is  appropriate  when  the  sheath  is 
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collisionless  and  the  latter  when  the  sheath  is  collision  dominated.  It 
would  be  of  interest  to  use  the  thermal  flux  relation,  plus  a  free  fall 
relationship  to  derive  a  load-line  type  equation.  The  equation  for  the  ion 
velocity  replacing  v^  =  &.E  becomes 

l. 

2  2 

vi  =  -(-SeV/r^  +  vs  )  ^22) 

where  vs  is  the  thermal  speed  at  the  outer  edge  of  the  sheath  and  m^  is  the 
mass  of  an  ion.  Combining  Equations  (5)  and  (22)  produces  a  free  molecular 

g 

load-line  equation. 

.  [2  Va  wi  (1  a«)]2 

3th  =  £s2  (23) 


2  1 
where  V  =  - 


nuv. 


■is 


2  k  To 


1  or  v2  .  -eVw 

=  (1  +  7^)  Jth  =  e  NsVs  is 

0 


,  h  Jt  Ok  \ 

'  1  e2 


the  thermal  current  density,  N  is  the  ion  number  density, 
is  the  ion  debye  length,  the  subscript  s  refers  to  conditions  at  the  outer 
edge  of  the  sheath,  k  is  the  Boltzman  factor,  and  Vw  is  the  applied  potential. 
In  deriving  Eqn.  (23)  the  same  boundary  conditions  on  E  and  V  as  before  have 
been  applied. 


When  the  only  source  of  current  for  the  sheath  is  the  thermal  flux, 
j  =  and  Eqn.  ^3 )  can  be  used  t'->  estimate  the  sheath  thickness.  Plots 
of  vs.  y  and  vs.  y  are  shown  on  fig.  22  at  z  =  270  Kft.  Since 

is  at  most  a  factor  of  3  greater  than  j  ,  neither  mechanism  is  sufficient 
to  supply  the  experimentally  measured  current.  But,  at  lower  altitudes, 
is  an  order  of  magnitude  greater  than  and  hence  in  rough  agreement  with 
the  experimentally  measured  current.  This  is  consistent  with  the  findings 
of  Hayes.''' 
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Higher  altitude  results  shown  on  Fig.  22  indicate  the  thermal  current 
to  be  no  more  than  a  factor  of  two  larger  than  the  diffusion  current  in  regions 
where  the  diffusive  flux  is  toward  the  probe.  Here  the  thermal  flux  is  cal¬ 


culated  using  the  simple  kinetic  theory  expression  for  the  thermal  velocity 


where  T  is  the  ion  temperature  =  neutral  temperature 
m.;  is  the  ion  mass 


when  Te  =  T.  Hence,  if  the  electrons  are  in  thermal  equilibrium  with  the 
ions,  one  cannot  get  a  sufficient  current  to  match  the  data  from  a  ther¬ 
mal  flux  without  large  fringing  effects. 


The  fact  that  the  calculations  of  Ref.  1  indicate  the  probe  inferred 
number  density  is  higher  than  the  predicted  N.  at  higher  altitudes  is  con- 
sistent  with  the  possibility  that  convection  is  important  at  high  altitudes. 
To  see  this  mathematically,  if  jgXp  is  the  experimental  current  density  and 
is  taken  as  constant  across  the  sheath; 


(26) 
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is  the  expression  for  HI  assuming  no  convection  but  just  a  thermal 
If,  in  addition  to  this,  a  convective  flux  is  assumed  to  transport 
into  the  sheath, 


^conv 


e  u 


(27) 


flux. 

ions 


is  the  convective  current  density. 


With  both  a  thermal  drift  and  a  convective  flux  transporting  ions  into 
the  sheath,  the  inferred  number  density  will  be 


exp 


i  e  vth+  eu  if 


(28) 


which  is  less  than  that  predicted  by  Eqn.  (26). 


It  should  be  noted  that  if  convection  is  important,  one  can  never  dis¬ 
cover  this  by  using  a  strictly  one-dimensional  method.  Instead  the  one¬ 
dimensional  approaches  will  require  the  sheath  to  becane  thinner  the  higher 
the  current  density  thus  rendering  convection  to  be  less  important.  If  con¬ 
vection  is  important  and  the  sheath  thickness  is  of  the  order  of  the  probe 
dimension  or  greater,  then  the  problem  is  strongly  two-dimensional.  Simple 
analyses  such  as  those  discussed  in  this  report  can  only  account  for  the  two- 
dimensional  effects  very  crudely.  To  analyze  the  problem  more  accurately 

would  require  numerical  solutions  of  the  relevant  equations  such  as  those 

7 

developed  by  Boyer  et  al.  These  massive  computer  codes  are  very  cumber¬ 
some  to  work  with  and  are  costly  because  of  the  large  amount  of  computer 
time  involved. 
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The  expression  used  for  the  ionic  mobility  in  the  load-line  equation 
was  a  curve  fit  of  either  Mason’s10  or  Dukowicz11  predictions  for  ionic 
mobility.  Mason  assumed  the  mobility  to  be  dominated  by  elastic  collisions 


and  resonant  charge  exchange  collisions  between  ions  and  neutrals.  Hie 
-fi  -12 

short  range  r“°  or  r  repulsive  potential  should  dominate  at  high  tem¬ 
peratures.  Mason  found  that  up  to  2000°K  and  Dukowicz  found  that  up  to  5000°K 
the  reduced  potential  Kq  was  in  the  neighborhood  of  4  cm  /v-sec.  Hence,  the 


potential  at  a  density  p  is 


K  =  K 


jr  - k  Up 


v-sec 


where  p  q  is  the  standard  density  of  air. 


Their  calculations  are  restricted  to  low  electric  fields  so 


E/p  «2  - XSiis 

cm-mm  Hg 


(30) 


If  the  above  condition  is  violated  then  the  ion  drift  velocity  v^  is  no 
longer  simply  proportional  to  the  electric  field. 


Calculations  at  270  Kft  using  the  results  shown  in  Fig.  23  indicate 
that  inequality  (30 )  is  not  satisfied  throughout  the  sheath  even  if  the 
sheath  fills  the  boundary  layer.  For  a  sheath  thickness  of  1  cm  at 
s  =  .37  ft,  E(0)/P  =  7*8  volts/cm-mm  Hg.  Similar  conclusions  are  reached 
for  other  altitudes  for  sheath  thicknesses  like  those  shown  in  Fig.  21 . 
Hence,  the  mobility- limited  drift  relation  is  not  valid  strictly  speak¬ 
ing,  But,  it  still  should  provide  a  first  approximation  to  the  ion  drift 
velocity  since  E/P  will  be  less  than  2  volts/cm-mm  Kg  over  much  of  the 
sheath  even  if  it  exceeds  this  value  at  the  wall. 
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IV.  CONCLUSIONS 


1.  Except  for  current  traces  associated  with  those  probes  located  on 
the  nose  cap  at  270  Kft,  the  experimental  results  were  hound  from  above  by 
the  total  convective  current  in  the  boundary  layer  and  from  below  by  the 
Q-l-D  or  diffusion  currents. 

2.  Since  the  experimental  currents  on  the  nose  region  at  270  Kft 
are  larger  than  the  total  currents  being  convected  through  the  boundary 
layer,  the  predicted  ion  number  densities  are  probably  too  low. 

3.  In  the  regions  where  convective  effects  are  important,  the  ion 
sheath  structure  is  too  complex  to  properly  describe  by  simple  theories. 

One  must  use  numerical  methods  to  adequately  describe  the  flow  field. 

If.  At  270  Kft,  on  the  conical  portion  of  the  vehicle,  it  appears  that 
the  inclusion  of  convective  effects  using  the  Q-l-D  theory  can  explain  the 
experimental  results. 

5.  At  no  altitude  between  220  Kft  and  270  Kft  on  the  nose,  did  either 
simple  convection  or  diffusion  theories  produce  results  consistent  with  the 
data.  But,  it  appears  that  an  enhanced  diffusion  current  due  to  local  charge 
depletion  in  the  boundary  layer  near  the  outer  edge  of  the  sheath  may  be 
sufficient  to  explain  the  experimental  results.  Ifore  effort  is  required  to 
clarify  this  point. 
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Appendix 


v 


Description  of  Load  Line  Program 

Card  # 

1  Title  80  columns 

2  Te(°K),  ue(gL),  Vv  (value),  FACRH0,  FACY,  NY,  NYS  (5F10. 0,215) 

(  1  if  density  table  is  gm/em^ 

FACRH0  =  J 

v .551  if  density  table  is  s lug/ft ^ 

C  1  if  y  table  is  in  cm 
FACY  -  / 

30.48  if  y  tabl  j  is  in  ft 
NY  =  number  of  table  entries  (30  typically) 

NYS  =  number  (odd) intervals  to  compute  sheath  properties  (i.e.,  boun¬ 
dary  layer  is  divided  into  NYS  strips) 

NY  Cards 

y’  f  >  e  =  ^  '*  Ne(-^-),  ~  ;  (5F  12.0) 

cm  ue 

YS,  NYPR  (F10.0,I5) 

=  ysheath  =  original  sheath  thickness  in  cm  (just  less  than  y,  ) 
NYPR  =  number  of  printout  station  through  sheath  (odd). 

Units  are  c.g.s. 


-23- 


1 


I 


t 


o  d 
o  o  o  o 
3*  0*  -I 
'  o  o  o  a 

►-►"KM 

oood 
o  o  o  a 

JIftP! 

o  o  o  o 

till 

o  o  o  o 
o  o  o  a 

•  •  •  * 


o  o  o  a  > 

«  •  •  •  Qmm 

•  ^  HfW 


uj  x  •  o  o **  I  oooo  •  —  -»  !g 

*  *  —  4S.  uj  o  .  O'-  •••«.  -*oc  it  t  a 

•  u.  •  (>**  o  ,  o  ac  ■  «  a  *  *  o 

©  ►  o  ►  •  —  GC  O  «-  ’  O  MMKM  >*  »  I  >  M 

♦  •  k  c>  a  cl  oj.to  o  o  a  o  a  >-  I,  K  o  o  o 

» (V  •  i  >  >  X©  V  V  •  O  O  c  «  •  •  •  Gmm  >*  J  c  K 

Ow  O  >  AUitfO  O  Z  A  HN  >  i  HN  ©  • 

nnoo  i  m  i  •  i  i©  i  «oo©ujhnmnu«w  ^ouj««  •  w 

«•*  D  in  —  —  —•«—•  1/1  cd  C*  ^  ■  con  X— .  —  —  —  OXK  in  Ort  D>>  N  m  J 

dx  j  j  j^o  >a  ui  —  —  «•  s  j  j  j  jo  n  •  v  <o  v  u  z  m  i/>  >  n  >  io  • 

ocac  h  •  •  .  .^ov  •oo-.-a***ujujujm  «-->■  <tovo-»nmv>-o  x  — 

.Ho'ctKB:N"<N  n  o  O  ■  >  K  K  —  —  VO  O  u  ■*'. 

x  *  -j  —  —  — - ■  a:  .j  _j2T~*-  —  h  w>w  j»»mi  ii  ii  w 

n  ■  n©7  •  <<  <«ja-  >  a  u.  onKtuatuou.u.u.ujvvvKxu.ou.00  —  —  >>o  «  x  u. 

*)  o  o  o'  m  ;  o  ?  j  m  o  .  I  o  |  |  ©I  I 

*«*  I  .  |N»  !  9  i  <r  I  ! 


^•c> 

>  >  X  o 

►  AUitfV) 

I  I  II  I  • 
—  —  «  —  i/i 


tL  >  [• 

>*  »  |  > 

I,  K  3  O 

K  -J  —•—■*»  © 

>  •  «-•  rg  o 

— 1  >-ow»-  • 

in  o  «n  z>  >  k  rg 

co>  u  z > s y 
«  o  >•  o  —  ►  «>>•© 

OO  ■  >*  K  K  —  O  o  u 


•  oi/iw^K^^^wuyyyoojtjc*- 


rg  in  m  ■  * 

•  <m  »  •*  •) 
>0  —  x 

w  iTi  tfv 
UJ  <  in  in 
K  X  -J  *-* 

—  OC  _» 

«n«on 
jr  a.  u  o  c 


c  r  •*Nfnift^No-OH(urunNoi»«rijrtg,inooc»»Nn^tfigjyo»r! 

'ru>miniMJiuii/Mnirin<o<og)o«<oor,NKNNyN(Oa3cocaco(Oti«}0'C>(>o,oiO,o,oo,o 
ejoococciooooociuooooooocooooaooooooooc.ooooo©*-* 
oc-cuooaoouooaaooooooaooooouoooooooocjooooooQ 

i/»  lAvivitrt 


n  n  n  ♦  in  g> 
a  ©  o  o  o  o 


z  r  r  r  r 

m  %n  %n  m  wi 


■  3  *  3  <  >  5* 

i  i  -  1  S 

*«  *»  ^  -  i 

,  ?  ,-i  t  -  l  ~ 


O  o 


o  o  o 


o 


o 


4 


u  ;  3 
.1 


i  •->  c  c  c\  r- 

'a  !  C3i 
in  ►  O  >  j 

I  ^  O  Q  H  II  ; 

•  r  «  n  r*  i 

■OC3-S  DD 

«  q  m  m  o  _<  ©  m 


!  * 

rg  31 

i  •/>  >  © 

,  •  •  o  — 

— •  1  hOOH 

*"l  ►*•-«*«*> 

*  ►  •  I  IM  <  ' 

-  N  Jt^TCN^' 
«.  •»  »  »w^(V(U' 
<:  -•  «\  c  c*  #  i.  ( 

>  1/7  <1  M  4/7  •  #  ►-  l 

i  «  -j  m  »r  -)  ii 

i  r:  —k  f\,  — .  - 

:  >  tn  if*  o  m  —  • 

:  i  n  r»  ■  u  H  • 

-.  n  _j  _j  m  < 

3  ->  ~  -J  ~  * 

.cgn  nc— 

r^UCtnilCii/U 


u;  { 

mm 

0 c 

o: 

1 

CL 

IX 

•  >  1 

UJ 

UJ 

i  Z  i 

•  +  i 

•> 

40  1 

> 

> 

I  Z  ! 

Z 

z 

*  <  1 

•» 

in 

4/> 

1  <  1 

z 

z 

I-  1 

< 

< 

UJ 

,  X 

o 

a 

1  »"  tt 

UJ 

u 

i  r )  m  n  m 

:  i  2 

:  cz  <r  =:  < 

*—  il  II 

•  LU  —  U.-  — 

-  o  •-  c*  ~ 


r  *1  «.  <4 

7  U»  4/7  W  4/7 


<  O 

4/7  CO 

•  C  U7 

1  —  in  m 

>-*  r*\  -O 

•  <0  •  < 

f*  f  s 

UJ  —  *s  O  < 

©  *-  c,  c  o  : 

*  -J  ui  c  c  o  . 

0s  *  O  o  o  • 

-«  •-  —  •  eg  * 

•  li  o  O  m 

Ci  *-  i»  H  N  v  > 
— ►  —  cc  o: 
vj  o  *-  — •  a  c.  < 

• C  •  *>  >  r 

•«*  —i  C.  (X  0  Q  , 

|t  w  —  II  II 

—•  t~.  «r  r  Mivt 
UCK.  I/I«JU< 


rg  o  O' 

x  —  »*  of 

-j  <  <  a 

*4/5  4/7  > 
•-  II  II  Ci 
II  •  <1 

—  ^  CN.  Cl 
u  w  C 
o  ~  © 

N  k  »  » 

.-•  c*  O'  rg 

r  <  «r  m 

/m; 


<  I  im  * 

in  _i  u  — 

*  -  < 

0x0  !—  40  — 

C  —  —  >k  o 

c  *s  ^  —  o 

.*1/7,1  k-  * 

^  ~ 

•JO  -  >11 

o  o  o  — 

-*  o  o  ► 

*  •  <t  Ui  •» 

o  *r  4/»  o  — 

—  i  i  •  — 

<  «“  "**  ““  UJ 

m  f\*  — *  -j 

#  I  J  ► »- 

c  _J  ^  m  — 

C!  *  * 

o  <o  h  <0  ui  «a  «■* 

•  —  2  4/7 

I  in  •  4/7  m  *«r  k- 

—  «.N  •  N-^O  C 

ii  ii  ii .  ii  »  if  ii  nj 

~  *-  O 

—  C  —  f*'  o  ►  — 

>  >  tt  >  •  O  O  UJ 

r  o  *-  o 

c  c  Ci  il  w  ■— 

iiiur  uj  ir 

c  {jccuumx 


im^lAiiOKCDC'O-^CMfnsrinOf 


)  i-  a.  n  o  **  n  i 

’  n  n  (*>  v*  c  xr  < 


Zz ^ zz zzz?zr rz z 

)l/74/74/>4/7i/)l/74/74/74/ll/74/ll/ll/7l/7 


zzzzzzrzzz: 

4/7  4/7  V*  4/1  1/  1/7074/74/74/11 


I  I 

i  ;  1 


o  c  QHPtnHpt-^HHk •rgfMfv*nirgc\iCMfMrM»N4r'<r'  h  r  n  g  r  n  m  i»i  v*  »r  xj-  «?  xf  xf  s*«fgg«'tf'iMAif» 

ooooooooooocroooc-oooooaoc  C7  o  o  c.  c.  oocooo  co  o  oor-  cooor.  o 


Z  Z  ZZZZZ'ZZZZZ 

ll/7l/l4/l4/l4/l4/74/l4/li/74/l4/l4/l 


o  o 


o 

o 


o 

< 

a 


|22 


J 


i  i 


*  •  > 

in  x 

,  «  «-4 

r>  k  • 

<N  Q 

;  uinn 

*  fttj 

c>: 

C  Zn 

»  CXO  ' 

X  O  *■«  t 
m  *v  * 
k  m  u> 
n  *d 
•«  x. 
m  —  3 

•  ui 
mn  »  ^ 

>**</>  x 

:  »uj  v  c* 

<  Jo  !  • 

(  _j  i_j  n  uj 

hzxc  . 

'DIM  j 

•  co  •  cm  , 

■  cx«>x  | 

!UN  I 

!  -J  •<«>«-«  , 

.  •>  •  Ui  < 
l  ©  <*l  ■  »- 
*H  juj 

«ujj:  o 

'<  •  2*“  —  i 

Ctult  ft 

X  X  x  X  » 

a  CM  O  -> 

■  »X  •  ft  x  ft 

UiWIOCQ  !«ft 

3  •  *  X  **  ft. 

•  ift  w  c:  u  • 

•u  «-«r  i  ,  «■» 

h  fn  |y  n  ^  w 

<  ■  x  x  < 

O  *unc?  »✓> 

c  >  II  y  N  •  ft  tf\ 

X  X  _i  .  • 

*N2*»ZI  »«4  >4 

X  rg  x  «0  x  x  (  NOiO 

ftJN  X  JM  ft-ft.  ft^  m> 

■ft  >0  w  fv,  ft*  ft.  ft  OJ  li, 

*r  ~  ^  c  —  o 

<L-<  HU  <  q  ■* 

1.^5  "I*-  X  >-  Ui 

K  ft*  C£  Q  (ft 

n«  oa^m 
l4,««Uu  CJ3UuOu>U. 

.  ni  m  i 

-4  cm  * 

O  O 


I  i 


WWW 


W</)VI 

i 


1  '  X 

Ui 

«  (h  O  *■«  N 
WW  o  w  Otf 

MM  «4Q 

c.  o  o  o  a  u. 

ZZZZZM 
ww w wwz 

ftft  ft*  ft.  ft.  ft.  O 

I  a 

< 


0  o  0  0  0 


0  0 


I  I 


M  ;  < 

.  »i  » 


0 


I  1 

:&i 

°>  ! 


i  • 


-23e- 


-23f- 


REFERENCES 


1.  Hayes,  D.  T. ,  "Electrostatic  Probe  Measurements  of  Flow  Field  Char¬ 
acteristics  of  a  Blunt  Body  Reentry  Vehicle,"  AIAA  Paper  $72-69!, 

June  1972. 

2.  Lew,  Henry,  G.,  "Shock  Layer  Ionization  at  High  Altitudes,"  (Final 
Report  on  Contract  F19- 628- 69- C-0LL2),  AFCRL-1 70-0702,  GE  7OSD782, 

The  General  Electric  Co.,  Valley  Forge,  Pa.  19101,  November  1970. 

3.  Bredfeldt,  et  al,  "Boundary- Layer  Ion  Density  Profiles  as  Measured 
by  Electrostatic  Probes,"  AIAA  J.  5,  January  1967*  PP*  91-98. 

!.  Su,  C.  G.,  "Compressible  Plasma  Flow  Over  a  Biased  Body,"  AIAA  J., 

3,  May  1965,  pp.  842-848. 

5.  Dukowicz,  J.,  "Theory  of  Convection- Conduction  Dominated  Electro¬ 
static  Probes :  Numerical  Solutions  of  the  Two-Dimensional  Flat  Plate 
Problem,"  CAL  Rep.  No.  RA-264l-Y-l,  (1969). 

6.  Johnson,  R.  A.  and  deBoer,  P.  C.  T.,  "Theory  of  Ion  Boundary  Layers," 
AIAA  J.,  10,  May  1972,  pp.  664-670. 

7.  Boyer,  et  al,  "Experimental  and  Numerical  Studies  of  Flush-Mounted 
ELectrostatic  Probes  in  Hypersonic  Ionized  Flows,"  AIAA  10th  Aero¬ 
space  Sciences  Meeting,  San  Diego,  Calif.,  January  1972. 

8.  Reilly,  J.,  "Probe  Data  Verification,"  in  "Research  in  Reentry  Physics/ 
Final  Technical  Progress  Report  (u),"  Doc.  No.  AERL  71-150  (December 
1970),  pp.  89-103,  Secret. 


-2  4- 


9.  deBoer,  P.  C.  T.  and  Johnson,  R.  A.,  "Theory  of  Flat-Plate  Ion-Density 
Protes,"  fhys.  Fluids  11,  909-911  (1968). 

10.  Mason,  E.  A.,  Estimated  Ion  Mobilities  for  Some  Air  Constituents," 
Planet  Space  Sci.,  1970,  18,  pp.  137- lMu 

11.  Dukowicz,  John,  "Mobility  of  N0+  Ions  in  Air,"  AIAA.  J,  8,  April  1970, 

pp.  827-828. 


-25- 


f/G,  7 


COM£>A/e/SOA/  S£7U/££A/  /tf£4SOG££>  SC)ep 
AA/O  P/?£D/Cr££>  CVgttA/r  £>£A/S/r/£S 


32- 


R  E  O  Jpjrp 
'SiTJES 


/Sr#S  A  CO^SA/r  D£A/S/ry’  J- 

*s  a  /reverses/  as  o/Sr*/vce  stou  a>*o3c  sSissci  y 


r/s.//  n/eoaer/c/tur  P*e0/creo  cazee/vr  &cA/s/ry  j 
*  fMcr/e/v  vs  £>/sr/wce  r*o/vr  />*oa£  svwce  y 


RG.  /S  7HE0*€r/C4CLr  P*GO/Cr££>  CO£p£A/r  DCMS/rf  J~ 

PS  A  W/YCVOM  ar&sr/wcefl&Al  M06£  SiMPce  y 


220  AC  FT 


o.  2/7  fr. 


AVCO, 


-9 


/r/<?-  20  &e/A/G//vG  careers 


